THE NAVIER-STOKES EQUATIONS FOR COMPRESSIBLE FLOWS. 
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y Ne: specific botal energy 
ya ME: fore 

F,G,H,U : vettor fuisce 

T: echian of coordinate trans formation 


Ri molecular heck conductivity 


Pr: turbulent Prandtl number 


М: eddy راصنا‎ 
A: bulk. масову Pent 
5,9, 1: transformed coordinates 
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1. Introduction . 


The twee dimenional , time-dependent , compressible Name „ Stokes’ equations ant 
Me cornerstone of Fluid dynamics. Applied to a control volume V eme by a contet 
عم سد‎ А, thue equations expun Aha conmenwation of maw, momentum and 
However Bren equations in thën integral form ane nob very usable, That is why, ымыр 
the transpre theorems , we transform. Hem into differential equations. Moreover these 
equation, ала exprimed dn a coordinate system. , generally the Cantewan coordinate system 
dn which the деа of there equation da much mou simpler. than in any ether coordine 
deme, and лол еп in a flux vector’ form which is the baic equation tn the 
numerical podus b be the Navier-Stokes’ equations - However , when studying, the 
flow struchins around a body such as ол cireratt , the Cantesan coordinates ane 
adequate in danibing the geometric contiqurations of the Girme and comequently ik 
and derive bhe. Arma етае отъ form of the Name, Stokes? equations, 


Кеййн ln the diusir aç Cox a a body such a» on aircraft , experiment 
haa shown Ahab ib 4 not always necessary Ло alle the complete thue dimensional, time depende 


compre " 
equations often give vow satisfactory nelik and enable to save commiderable computations 
atime Depending on the zone of study у at ù н бо apply some or other approximate 


an unsteady phenomenon arda у the advantage of the zonal approach over the 
fll equation зорго te uncertain. This wilh nob be discussed in the present paper 
buk let us remember that the zonal method le a vast subject ef research. In the 
rende betwen the complete Nowin. Stones’ equation and the beundaw foyer equabiens, 
the “approximate Naver, Stoke’ equations met wad for compremitle flous аге; 

= The “thin. * equations 

— The “panabetized” Маме. Stoku’ equations 

— The “emicat” Nave. Stoku’ equations 

— The Маси Shock. Layer equations 


There are many numerical solutions to the Сотырлам е, Naver. Stoke’ Equations. 
Howaven in the peni paper , only Ahe finite. difference proudu ane addremed. 


2 
Beam and Warming 's approximate factnizabion achame у the Briley and MacDonald's 
procedure and seterad hybrid methode, In this puent stent , the scheme of an ici 
methods and the scheme of an implicit method are given, the Explicit MacCormack.’s method 


The principal aim of this paper Да to introduce the neaden to Ahe numerical 
Ашота of the compumible Navier, Stekes’ equations ; 
= Маб à the general form of the equations aec wish to sele? 
— How and for which Н ima these equations can be simplified ? 
— Which finite difference schema ane applicable to there equahlon ? 


2. Derivation of the compressible Navier. Stokes’ Equations. 
=з = ca alkis Se e ыы аан Д 


26 In &he Eulerian formulation у the conservative equations in integral foun for 
таз, momentum and energy with I Жо a control volume V atationony Ал Ahe 
dritioh frame and enclored by the control surface A ane: 


e |c) EXE dV + TEE dA = о er 3 iil ev + flea dA =O 
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ee m ыы enî 
e The second equation i the Newton's second four of motion : 


eje eee 


The умло of comowation of energy (equation (3)) states Hak: 


fole of increase flox af energy work per unit of time * heok flux 
f total energy inv = entering V accross A)” e the pressure and viscous forces) (entering V accross A 


a iege ы iS Generated } 
af the body forces inside the volume V 


Now, since Ае control volume is stationary , he relative order of bime derivabin 
and vere integral. may ba exchanged у ond by means ef the divergence theorem we 
can derive. the corres pon ding differential equation ; 


e Be + Vlesja o or + 2 eu) =o 


* (5) Fes + 7, (el + p§-5) - ef =9 
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In thee equations , the viscous stu kerim is given оз: 
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let w change now our notation in the following way: 


5 — (vw) > (B, Z) — (xyz) 
And Ды u introduce equations (7) and (t) in khe equations (4), (5,), (52), ( 
and (6). We get a system of five equations where ther ane ten unknowns 
(4%), e, p, T, e, К, 9, A). Therefore additional equations must be included. 


For mæt asnodynamics applications we adi khe following equations : 


# İdil) The Sutherland’s viscosity Дош; 


32 


T+ Ca 


e |002) The Prandt number à assumed to be constont : Pa SEP 


? e | (13) The Stoke's hypothesis (zer bulk Vibeotiby) is assumed to be даалай : 


3A+2p 20 


Thee five addicimal equations complete tha deroription of the system of equabons 


POP 


We have 4= - VT and ce 
Cp 
is Th maa ў@ шс M 
е 4 a тае 
o Тл Ае equation of energy expressed in terms of the partial derivatives of e 
E OE 


At this step of our discusion , we now he bod f 2^ the 
external heat addition û. Я 


We must alo comidos separately the case sohen the flow is Laminar ond the 
case when the flow ù turbulent. Е. 


heat s 

Then | Bo: inte approach.» We can change the Naver- Stoker Equabiow 
be a тобед form of the Reynolds. aw equations by replacing, Але соге 
pone она dü о а E mein скн сы; 


Mr © the eddy усов Uy - ± (7.5) § + М (VF +907) = ЕТ 


kp e turbulent thermal conduchvied ; oR, VT = 


Y turb 
_ Then me can dine te tunbulenk Brandt number by ele 
Р 


And В. У.Р тё 
Consequenkey д, = y ve FE ED 


+ In Are equation of energy expused in terms cf Айе parkal derivabve op 
ё , we сол лао ğ ey OBS Ву 
D 


- The Nowier. Stokes’ Equations for Laminar flow and Reynolds. ow Nav 

ол. + overagad - 

Stola? Equabions anc mest jrequant®y numiten in vector form болул. 
Coordinakes , This form à: 
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(ee * p) СЕБЕ — (ux + V Tay + ит) 
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These equations ane denived with aspect Ao Khe Cantertan coordinates. But for 


GD cing реу System of саа рот ERE E 
the qtometuic confiquration . First because an extensive interpolation procedure 
E cape ions, ааа ao mame stieg oA 
© systematic 7 id spacing . Thus a У dinate ن‎ b 
A сеа шы Te р e? 

== (8,0,5) y= ¥(3,9)3) z=2(8,9,7) 


Then equables (14) becomes: 
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where J= enn). Mu ü the Jacobian of the coordinate transfor 
(x,y,z) 25,9 3) 


and co LV ek transpoie of the abe of the maki Playz) 
2(8, 9,3) 25.92) 


Tua Nx = == 75% j р ур - уду ) f. Ys 75 =- Yo Zg 
SS = тұзу - туту ). уре Eea зү лүш ар: 
3/7 = 71-77, ) Pals = =z Ys -75 Y7 j lasti- 
(io gute Сы eh, Bos анына нр адо: 
adding, and subtracting Hike terms , We have: 
(2) скен » + Sov e «(Beri Bie) 
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The tat thin terms ane all equal bo zeo and tan be dropped. And 
Sua ر‎ the squobons in the new System of Coordinabes take Aha abong conserration, 


PET $3) & aue By MES Ais e Gp, + Hy, )/3] 


+ 2 [0 + ау + Н5.)/7] =O ur SN 


Thus we derived the strong comevation fom for the compumible (lamin ov T 
б кален че асое ат пасри eles i 
= zero become we ragltcted Bu body foru £ and а 
we had not neglected icc they would haue appeared in the night hand side of 
thue equotions . 


Note: » The partial deme tems appearing in U,F,G ard Н ane Же be Eranatormad 
heo. For example Ahe viscous strum sensor sould be tranformad to + 


Tia Que = Dus Qu; 2ч; 
= у um ur келу їр 
: Gps n, MTS lege 


e Ij mu girmeli, we had omumed a Coordinate transformabion. in the form : 
=el, 5, 3,%) y:y (G9, 3e) z2z($,9,2,5) 
Then Ahe strong Сотые от, foun sould Фаме. taken dhe form: 


20/3) + (03. + F5 + АБ, = ып + so Fy, + Gp, + Hy, )/T] 
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3. Approximate Navin - Stokes’ Equatims . 8 


А Tk 0 vey difficult to solve the compute Nawier. Stokes’ Equations , buk 


) fortunately , {оъ many studies of fous around a body where she boundary. Фауел 
equations ane поб anplicable , it ip however possible to solve a reduced set of equations 


that bekuan the compite Novie. Siu! equations and the boundary Layer equation 
an terme of complexity. most frequency used approximate Nawier- Stokes’ Equations are: 
22 = The “thin. Y fexto: Stolus! equations 

=W t ud” Namen. Stoku’ equabons 


s The “conical” ~ Naven- Stokes’ 'equabioms as „у ha 


Dd RP )* £v» 


(20) | у momentum: We Rew) Rare f IR (evw)=o 
% momentum: BE „е es a м5) +2, (p+ ew) =o 


EE 2 2 
energy : Аааа и): ey gey p eme serin Lİ 


+2 (gew + pw) zo 


lak ш note that vàcou, terms هه‎ 


well in 
fux pine geyve. Thexfou we ninin only the heot in Өш. у. dinski i-e 26 _ 


= 7=#ў(®›у,®) 
3-3 (= ук) 


Rudman and Rubin (1968) wena олом the firsts до dee the PNS equations 
Their y 


(0) nondimenwonatize Ahe flow variables with nupect bo the freestream value and 
the space coordinates with nupeck to chanacteudtic лаһ (Lin the x-ditechon, Б 
an the y. direction) : 


ьа ож Va у= v? Me ps P“ Pa e= e*e. тетет PY vs 
Race” LL. ysy* ad  Sz85*L 


(Ù) омште series expansions of the nondimensionalize flow variable: with uapeck to faak 
flew variables mondimennionalized wiih respect to Local ruaference conditions : 


Р = Pree (I + ере +---) 

Career (ee eee) 

T's Tee (Tot + етю...) 
Е, 


Ps ра СМ + є tes) 


e The zeroth - order solution d obtained from the first term in the series expansions 
while the first-order solution ù obtained from both the fist and second term in the senis 


Furthermore dn the relatively thin üzüm where the flow à disturbed , the gradients 
normal to the surface are much greater thon the gradients parallel to the surface and 
S^ can be assumed to le small. 


(ili) substitute the expanniona into the two. dimensonal, steady, Namer. Stokes equa - 
tions and do a magnibide analysis . 
The nk i» that : 
— The zeroth-order equation are denved by neglecting berms of order 6*9*, 
Ma Tre / Ma у end €. According te Rudman and Rubin çin order for (5%)? to 
be vay small (ie. 0.08) the zeroth-order equahions are not valid Upstream. af 
the point ak which: 
= Ta /2 Mİ Re “2 
легна where Ж e] ("la Raze) 


The zeroth - order equations are alo applicable only when Ma 75 


continuity : Pet + Ber „О 
pemenn (ye reu + SS (WR) UEUS- E202) 
energy : enc, dL 4 е, 31 = -р (25 КЗЫ ыу + rg + к 


— The first-order equabions ала applicable when Me z2 : 


and 5:5, БЫ ч 


The zeroth - onder equabiony become : 
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Continuity: Peu + Феу + dev O 
эх dy ez 
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x- momentum (go сы +еу a, On Ed = & 
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ya momentum; е“ 


m 2 
$$ (2E + Bt) + 


dw 2 2. S 
Z. momentum: Bis rev ers SAS 4 (м 55) ا‎ 80) + & (vos 


еч c, т чс us 25). 2. > (кэт) 
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Lek ш note that in the PNS equations denived by Rudman and Rubin , there 
à no streamwise nasus gradient term and consequently the equation Gin 
an a sin, “panabolic” manner 4n the boundany. Layer sugiom ion. 

Another form of the ENS адыны: dad e Gi TRES. PUE eal 
benm -9p/9x , and probably the most common form of tha PNS equations is 
obtained by neglecting the streamwise Views derivative tive terms у Седа heat _ 
flax terms , in O(I) compare to the normal and олмеш. [ein ce S, ù 
terme in O(Re3) . The equakion dn æ Cantestan coordinata: 5 System (x b the 
streamwise direction.) ane: 


continuity : same as equation (21) 


Sec meleke Come Sete be ein eC > (95) + & (8 52) 


Эх. 
res «фе ven “BF 880) özle) i Ug) Ê 
9 2.( 2 p2 
2. momentun : eu BY «ev o m зар (vt) + Bw) +203 22)- EE 3 


energy : same as equobon (21) 


Th order to apy the above equations to any more complex body sufoce , we 
must now consider the PNS equation in thein strong conerwation form jie. Ал. 
bam cf a generalized coordinate system. Recalling that the flow à steady, we 
can wite: Н 


AIR- 9I + (&;- 6E, + (Hi - 81/2] öper. * Саа, ei 
*alh-5-* (&-6,)1y + -H92,/7] 5 


where the flux vectors F,G and H have ben decompouc in thein iwacid and viscous componer 


coo/ece 
| 


Then neglecting khe streamwise vacou denivakive terms у the strong conservation |2 


FB, + GLB, + HE (ж streamwise direction 
Y z ) 


" 


((F =F), + (&;-67)p, + (Hz - HY) 92) 


;- КЗ + (8@;- 8), + (8; 82)2,) 


and Fy, @, and Н), are the viscous components of the {lux vectors F, G and Н 

dn which the terms containing partial derivabives with лере to E have ben omitted, ` 
For example Tans Tay, Tez and Gy pear in F, ола M. 
фла 2 p [2{5x + 9 Uy + £y) - (5y Ve + Dy My + Зуу )- (E Wt Dp + 3z“ ([ 


Tay = (By uy + Dy Uy + By Uy + фу аур * Zag) 


Tua (bz vg +9, „+ oz + BM, + Je Yo + чү) 
а = RE Ty +9, T) + 327%) } 


e These term an changed in Th, ET Tz and qj in F} where; 


Жа $r I E чу)- (Dy у умр) - Oo + 3, ;)] 


ту = PO. + 3 9z +9 و‎ + vo) 


Та = (9.9 + 3,674 Da Wo * 5. м) 


== -kO + 2, T7) 
Р if We have : 
eu ev ew 
2 езу euw 
p pos G; = [оир Н; = [суу 
бом uw wisp 
e+p) (ee+P)y (ee p] 
o o o 
Txx тку xz 
e Gye] wy Hye! tye 
Tez tyz Taz 
UT.k + xy *WIxz-4& “Tay AV Tyy + уту, -4y "i tYTy pt WT. "T 
rN This tema couewation foun in terms ap We «лудо and viscous pan of the diua vector, 
x tb employed да numerical prowdunes for seng the PNS equations. 


Th many studies , the “ thin. Loyer approximation.” can abso be applied to the PNS 


E equabioms . The nnutting equations ane tre steady form of the thin - ауел. Nawier- Stokes’ 
a ime _ á eon ob na. 2, 


3.3. The “conical? Navier. Stokes? Equations . 


A сото? gow à duch that all slow properties along the radial direction 
an constant and еласа a conical flow problem is a two. dimensional problem. 
Buz 4b ù dvi that the above property implies that the flow be inwiscid and thus 
the соле of comical flow is strictly valid омду for Асма inae, However 
experiments proved Phat the viscous portion of a flowfield sar strongly dominated by the 
cuter Anviscid conical flour. Comequently ر‎ according to Andersen (1973), one can quickly 
timate the heal transfer and shin frickin for such flowfields by solving the 
unsteady Nawier. Stokes’ Equations in a time- dependent pendunt on a sphere with 

Молу flows around conical or pointed new. d a have been computed 
eee hae sgn NEN UN уаш the лаһ; provide балла 

+ 4e- PNS calullations , 

The “conical Мамет. Stoha,’ equabions ans derived from Айе complote nondimen. ` 


“dionalizad Nawier- Stokes ecquabions "s 
РОР те (BES у 209 0. 9U'* T S 
2t axs Dy” az” 


where the ux vectors U*,F*,&* and Н" have exactly 4he same forms as those 
ben in equakons (15) to (18) by replacing Ае flow wuabls by their nondimenwona. 
лад countrpants which ane: 

ТИ ГЫЧА wh = мМ rs FI Ma €" Eleu p“= ve/e, 
тту e = е/2 j &*sz/L ‚ y*=y/L z*iz/L Me eV, /L 
Then defining : A 

Marta (x024 yM e uti) و‎ pee ў TE j Asie phe n у та E" us 
«we hoe the following, arma comsenvabion form: 

a (st M / 

e (3 YE (Ec хране) SUR carte ач) 2. (Ss Cate H*)) =o 


Эх. Эс. Әм. 


p ; Then aneurin, SET, 268*,.24*.6 and compubing a actubion or a ephe 
i whore nondimenaonaL nadis Да ле а л/а D xol we obtain the “conical” Name. 


(24) 
İn thue equabions , the partial derivative appeaning 4n the vibcours terme 
q F', G* and H* must be transformed Алла) : 
. A а 2. 2 h axa 2.3. А. 
axe = - PAGE Vet 2 5» = > 2 2z* Aa 
э» lt us note that tu “ conical” Nowier. Stoka equabions depend on the porikon. 


nal whee they ала computed though the Reynolds number Re, « o Val which 
appas in the vaou benm, For iwùcid solutions у thae equabiow 
ora independen of n and thus Ачу comical. 

seals 


+ atii e 
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The “cemal” Navin. Stou’ Equations are a Very approximate form of the 
compie Nowier Stoked Equations and comequintly an to be owsided for flow problems 
= wuwu a high degree af accuracy ù required. They take a very simple form and 


) can be solved with a time - dependent algorithm . 


5.4. The ¥acou Shock - Layer Equations « 


The Viscous Shock . layen Equations show many adwantages : 

— They are simpler. khan the PNS equations 

= Sine Ahay rumoin fypenbolic.panabobie in both He streamwite and 
cots flour dlinections , they can be solved by a marching pridu 4n. both 
Флит le the PNS equations mhich ant hyperbolic. panabelic only in the 
dtreamwiee dinechion mush be solved simultaneously over the елЫл. crouplour 

1 
NOU. терра MED ME СС 
PNS equations ала nek applicable since the flow ù subsonic in tha ragion. 
Conaequenkly they provide a stonbing aclubion, for a subsequent PNS resolubion . 


However , thein major disadvantage ts that ae eee 
flow fields m. emp s MART. ec = : endi a iin 


The Viscous Bilek Mops Equakions ana: 


de [(r* + 9% cos à)" e" ^] + öy [1+ K*p*)(r* *9" cas d)" e* v*]- o 


Continvity : 


E momentum: о ki Sigil o, ! Эь* 
I-K*o* agt . 27" \+ ке рн lek ЕТУ 
2 
є 
“kiye *)2 (Ma ay ap ОТ? *)2 (+ pM cas d) x*] 
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where ee ке ця ) 


Tak pF 
J momentum: e Шр i + ve Ov™ „Ко ыыы; Эр р 
\+ K"g* а 99* ї+к*о* 


energy? e'(Lut E et) ot - vee , (te)? 
ERE Oe ag% k 


Ca $5 +7°) +° ”ر2‎ ed 25) 


nnen 
C+ K¥ 5") (79 + пв)" Br Op" 


There equations ham ban expamed up to second order in € where 


1/2 
) and Мер ü the cæffiuen of viscasity evaluated at Trep a Vo 
Fu 


m=O corresponds b a 2.D intrinsic coordinate system and msl t» an 
axisymmetric. body intrinsic coordinate System: 


eoo. 4 


2.0 or exisymmetrie body 
intrinsic coordinates (1,9, 9) 


ЖЫЛ ГУ 9* = 9 / nose j 

ож. v/ Va j vw V/V j e" 

Tea T/T 5 v Pl Pree 
pt = K* e" (u9* 
СУМ | + қ 0% 


"S 5. Conclusion: summary оф the. results. al 


AFTERBOOY AND NEA 


WARE 
Necenity to have an appropriate turbulence 
model for describing the abrupt flow 

Structure chenge and its impact on 


WING « FUSELAGE 


FORESOOY 
Wide usage of the TLNS and 
PNS tmaotiom solved by the 
Finite difference methods mentio 
in Ch. 4 The ful N-S. equations 
were also applied fora comparative 


PROPULSION SYST 
“AND 
COMBUSTION 


by strong invistid -viatou. 
аа with a multi 
si waves System. 

А turbuionce model ما‎ айо 


The PNS equations have 
ben uted for this problem 
For subsorue problem , tha. 
Фу consists in 
implementin, brundary 
Conditions ; this diff 
{5 however compensated 
by the much belke 
undeztandied cf kurbuler 


AS boundary layer 
ناما‎ yer (subsonic) 


BLUNT NOSE 
The resolution af the VSL 

‘equations providis a starting 
solution fore subsequent PNS 


However CNS equation ore 
te be used only if a high 


degree of accuracy i not 
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Newtonian Fluids: Stress is a linear function of the instantaneous 
velocity gradient. 
Non-Newtonian Fluids: Stress is a non-linear function of the 


hystory of the deformation gradient (fluids with "memory": honey, 


heavy oils, molten Plastic and metals, jellies). 

Entropy: Measures the excess amount of heat generated in a process 

E beyond that whích can be recovered as work done by adiabatic 

expansion: EH 

-- slow compression heats the gas, while slow expansion to the 
original pressure brings it back to the original temperature 
(all the stored energy is retrieved). 

— fast compression produces compression shock waves that result 
in excess heating and produce the entropy. барбай to the 
original pressure leaves the gas at a higher temperature than 
originally present (not all of the input energy has been re- 
trieved). 

-- a fluid with high velocity fluctuations has more kinetic 
energy than one with the low velocity fluctuations. Vis- 
cosity acts in such a way as to decrease the frequency of 
velocity fluctuations thus decreasing the kinetic energy 
which then converts into heat (or internal energy). 

eL in the case of a compression shock the conservation laws are 
satisfied (mass, momentum, energy). If the shock is very 
thin, then the compression length and time scales are neg- 
ligibly small. Hence, there will be no velocity fluctuations 


after the shock, and this would mean that kinetic energy of 


m д ——— : ——M——: 


the flow decreased, internal energy increased, thus, entropy 


increased across the compression shock. 


NOTE: The system is insulated! “Therefore, if 4 = 0 (adiabatic 


ant eet 


walls of the system), the work Ф is not really dissipated and 


lost, but is retained (as an internal energy) by the system. 


C eva PA ЖАД 


a arbitrary base 


2 vectors 
ie) 
ES 
du - 
02 03 
ай ате es ах, ах, dX, 
аук PX) Boz By, “01 “oz оз 
lor, эх, эх, * 2 
Ap оза 
9X, ax, 3X, Я 
dV = det ах, ах, dX 
рео аиа. 01 402 03 


dV = J dV 
o 


Note also that 


ЭХ i ex 


is d аху OX, әх 
dt “ijk dt X 


G x E 
aX, 9X, OX, әх 


= є [stata at 2nd term + 3rd tera | 
ijkl ЭХ, 3X , ЭХ Fok 


But 


€ т mee 
ijk Эх р 9X , Xok 6823 


Hence, 


i 0 ax. ax : ian 
dJ 2 3 p Фф 
с С 35 ax coe Tes oy 1245 xí e 
d. * 
at ká Ñ u 


Let F be any arbitrary scalar function. Then 
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dt dt 
v 
o o 
dF dJ 
| дє ЗҮ, | Beas 
v 
o мкад 


Finally, the First Reynolds Transport Theorem: 


а dF ә ә 
$ | r av > | СЕ + Fea 
vy me. 


Note that it can be further expanded as 


dt ot D at 
у 
2. 2 а > 
= [tes (u - VF +F $-3)]av 
8 v 
Then 

a - (22+ $-otyjav - [ 3E Е 
з dt | F dv je + Ve (Fu) ]dv | ot dV + | F(uen)dS 
3 v ۴ M av 
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CONSERVATION OF MASS (CONTINUITY EQUATION) 


| av, = | pdv = | p зау, 
x у 


у 
o 
Then 
2 Lagrangian Form 
| (р -pJ) dV = 0 р =pJ of Continuity 
o o gres o 
v Equation 


o 


Otherwise, use Reynolds First Transport Theorem and, instead of the 
arbitrary function F , substitute fluid density р = (Xt) 


Then 


DAVIS fjg + p¥-ujav = ie + È- (où) Jav = 0 


у \ 


s 
په‎ 


Hence, the two most common forms of the mass conservation equation are: 


non-conservative form _ ё + pveu =0 o 
Эр ў. (оў 
conservative (divergence-free) form sc *(pu) = 0 


Note: If there are any sources or sinks of the mass inside the flow 


field, the right-hand side of these equations should account for them. 


Let È be an arbitrary function. Then 
Е d $ 4,2 
dt | FG dv | ёа А 
у M 


Now, use the Reynolds First Transport Theorem and get 


ss 


2 
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Гео + FG (Wu) Jav 


у 
+ dF © + > + 
= fe ЯЕ. рс, eösüylav 
E dt 
у 
тав MIC 
= fe (SE reu] + r Dav 
М = 0 (Reynolds First Transport Theorem) 
a Ж ac Reynolds Second Transport Theorem 
dt | FG dV = | Р de dv 


F =p and é = a where а is the fluid velocity vector. Then 
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CONSERVATION OF LINEAR MOMENTUM x 


Types of forces: 1) body forces: ob 
ii) surface forces: 


pressure (normal) * 
sheer (tangential) 


k a 
7 stress tensor ; T = I'n 


4 | sü av - f o av + {та 
у у эү 
Rate of Change of Body Forces Force of the Traction 
Linear Momentum (Gravity, Electro- Vector on V 
magnetic, Coriolis) (Pressure and Viscous 
Forces) 


Use the Reynolds Second Theorem on the first term and the divergence 


theorem on the last term and get 


> 
du * 
[^S f o av + та 
у у M 
or E 
di 
iu . 
ааа: ЖУ 
Since Е 
ТЕН сызгыра en; 
Hence, 


d 
p E * p(u-V)u * ўр - ob - -o= 0 Momentum Equation 
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CONSERVATION OF ENERGY 

First Law of Thermodynamics: Ae = 6q - êw 
ög = specific heat added to the system 
6w 7 specific work done on the system 


New Variables 


e 7 internal energy per unit mass 
a 22 
п = = kinetic energy per unit mass 
q = energy source per unit mass 
4 = heat transport vector [t1 =- k 97) 
Fourier's Law 
27 


$ 7 body force per unit mass 


22 
= | pi} + elav = | pb-u dV + Fus + [oiv = finas 
V V av V av 


Rate of Change of Work of Work of Rate of Heat Transported 
Kinetic plus Body Traction Heat Across the 
Internal Energy Forces Vector Generated Boundary of V 
(Potential Energy Internally into V by 
Neglected) Conduction 

AE -6% 8Q 


To transform this integral equation into a differential equation, 


first note that 


sv 


& - ads = fe Дау 


av у 


14, 
(298 *T) + uds s= m dv 
T * udS = ÓO(f-I) - udS = " Tij ч; d 3, ij ay) 
av av © av у 


[ON + T1۷ = İs * Сует) +D: TJav 
v X 


ðu 
Since 7,424 = Tij By = Tis (D, 44255) = "j =D: I = stress power 


ae 0 (I is symmetric tensoryad 2 ú ал aed ipee fete tate), 


Also, use Reynolds' Second Transport theorem to get 
+2 E 
d u + du, de 
dc | 015 + eldv | elu a arly 
¥ у 


Finally, balance of energy can be written as (if â = pq) 


dt 
у 2% у у 


x 5 
aco - o8 mave f odeav a fips gs û ê - dav 


= 0 (continuity of linear momentum) 


where 


Hence, p HE р: r d 9-4 Energy Equation 


The following variables are assumed to fully characterize the 

thermodynamic state of a fluid particle. 
P = density 
T = absolute temperature 
s = specific entropy 
e = specific internal energy 

© These variables are related by caloric equation of state 
e = е(р,ѕ) 


T = Т(0,5) 
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1 


16. 


The Ist Law of Thermodynamics is given by 


2 2 2 
| - [ow = [6 те, - е 
1 T 1 
A where q = heat transfer to a system (per unit mass) 


Ы м 7 work done by the system (per unit mass) 
e 7 energy of the system - kinetic energy 
* potential energy * all remaining types of energy 
(per unit mass) : 
The symbol 6 denotes that the quantity is a path fection and, therefore, 
6 is not a differential. 

Gibbs relation is derived by first assuming we have an inviscid 
fluid undergoing a reversible process. The inviscid fluid assumption 
then stipulates that all work done by the system can be done only by 
pressure forces (normal stress). ch addition, we neglect the body 


forces. Thus, 
E 1 
öw = pdv = P а) where v = ry = specific volume 


Next, the assumption of a reversible process implies that the entropy 


change is given by: 


Таз = ög 

n Thus, for a reversible process of an inviscid fluid, we have that 
1 1 2 

de, = Tds - Ра) = de + du) 


where ec specific total energy ы 


specific internal еперу 


e 
" 


fluid velocity magnitude. 


Since all quantitites are now in terms of exact differentials, we can 


"divide" by dt to get 


E a 


Gibbs relation was derived for an inviscid fluid (with no body forces) 
undergoing a reversible process. However, the equation contains only 
state properties of the fluid. That is, there are no path dependent 

variables. Thus, the equation can be applied to determine the change 


between two thermodynamic states, regardless of whether the intermedi- 


ate process is reversible or irreversible. 


Therefore, Gibbs relation is valid for all possible changes of 
state produced by either a reversible or irreversible process as long 
as the change of state is calculated by integrating Gibbs relation 


along a reversible path. 


Here we assumed that there was no other form of work being added to 


or removed from the system, such as shaft work, and work done on the 
fluid by, say, electromagnetic forces. If other types of work are 


present, Gibbs relation is modified to read 


de 
ds d l t 
er ene ar меча E E 


where the бе, › 2 5 1 5 п are contributions to the total work due 
to other processes than pressure forces. 


Assume further that each fluid particle satisfies Gibbs relation 


H 
de тз D EE puree 
dt de dt 22 t*Tar 
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Gibbs relation and the caloric equation of state imply that 


de _ Зе do , Зе ds . СЕЕ 
dt 9p dt ds dt os as 

LOT 27 я M үрә 
ат = 35 dp + 3 ds Ж. 30 7) 


Note that from conservation of mass 


> us uod Ll P PET. 
d teu 0 TN REC EE gum 
Then 
T 
a} 
Mog oe R с жае ан ee Dye, 
at 2dt р p D*$5PD:1 
p 
Hence, Gibbs relation becomes 
l 
a 
de. es -ri Pp, 
de Du Р. тає" аваси 


"Combine now Gibbs relation and the balance of energy 


de - ds } А + 
aE PIT a SP DE а kö li 
Then 
Pitan AC TEZ 
ora 2: Pete 9 


Entropy Generation Equation 


VISCOUS DISSIPATION FUNCTION @ : 


6-р: (r+ 1] 


@ Recall that 


D-ireres =- pl +S | where | tr S$ =0 


„ е 
Tensor § is called the stress deviator or extra stress (s consists of 


the off-diagonal terms of Т). 
Note: р = р(р,ѕ) = "thermodynamic" pressure 
P=- i tr I = "hydrostatic" or "mechanical" pressure. 
For the fluid in motion viscous dissipation function ¢ is 


=D: [(p-p)I + S] 


For fluid at rest: 


IMPORTANT: Viscous dissipation function Ф represents conversion of 


mechanical work to heat. 


Hence, entropy production equation can be written as 


Viscosity Heat Heat Transfer 
Source 


Dissipation Principle: The total entropy Sg of an adiabatically 


enclosed system cannot decrease. 


5 _ (50 + 
Ls ve * 45, prev 
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> > 
р СИЕ V^ Mİ de 8Q = dE + би 
dt G T 
у av First Law of 
Rate at Which Rate at Which Rate at Which Thermodynamics 
System Gains Entropy is Entropy Crosses 
Entropy Generated the System 
Internally Boundary 


Use divergence theorem and dissipation principle EA 2 0) 


Since V is arbitrary 


ds 0 lee 
o. MÜ ibi. 20 
But 
= - 21094 + 9.09) -itam + ف‎ 
Hence, 
ый CLAUSIUS-DUHEIM INEQUALITY 
7 20 (holds for all processes) 


MORE ON DISSIPATION FUNCTION: 


I ETE 


Stress: 3i 


(general expression) 1 - PI * 


‘e 


Rate of Strain: D = Dji D22 023 


2+ 
tr D= Dii + P22 + D33 Veu 
ðu ðu 
B D Ee Econ 
1j' 2 x BE 
j i 
Stress Deviator: S = AI tr D+ 2u D 
Since we know that T is symmetric trs=0 


Then 
ЗА Vea + Û =0 7. (ЗА + 20) (0-0) = 0 


Hence, if [224 # 0 (compressible fluid) 


А = = 


we 
B 


and we get 


Stress Deviator 
or 


s--2u1@-% +2 p 
Extra Stress 


3 


and 


Stress Tensor for Ыы 
Compressible Newtonian 


1 = ЕВЕ ов = 0-011 + 2u D 
Fluid in Motion 
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where р = p(p,s) (thermodynamic pressure). 
Hence, the difference between thermodynamic pressure, p , and the 
mechanical (hydraulic) pressure, p , is 


=, E 
P- P= vg (Vu) 


ug is "bulk" viscosity which, therefore, exists only when the 


compressible fluid is in motion. 
VISCOSITY COEFFICIENTS 

u = shear viscosity 

А = "secondary" or "dilatational" viscosity 

My = "bulk viscosity" (sometimes incorrectly called "secondary"). 
Bulk viscosity is related to the relaxation time (adjusting time) of 
energy transfer from translational to vibrational modes. Hence, it is 


more pronounced at higher frequencies. 


2 
NOTE: A = (uy - à) 
-- "Secondary" viscosity never multiplies off-diagonal terms of 
the strain tensor р. 


-- "Shear" viscosity multiplies off-diagonal terms of D. 


STOKES HYPOTHESIS 
He assumed (and later expressed a little faith in this assumption) 


that thermodynamic pressure 

de 
а ` 
© 


p(p,s) = 


is equal to the mechanical (mean hydraulic) pressure p = h tr I , that 


“is, he assumed that 


22. 


23% 


2 
lytA +0 


This is reasonably correct for monoatomic gases (neon, helium, xenon, 
krypton, argon), but is completely incorrect for most polyatomic gases 


and especially for bubbly liquids and foams. 


Consider incompressible flows when U = tr D = 0. Hence, shear 
viscosity, y must be non-negative in order to produce non-negative 
viscous dissipation function © . This implies that A 2 0 (always!). If 
this is true, then Stokes' hypothesis Gg = 0) suggests that 


А 5 0 always, that is, A= - 


we 
= 


Since ug (bulk viscosity) is a measure of a deviation from the Stokes 
relation 2 

W23v9*2 , 
It follows that ug 2 0 always if we have a purely dilatational motion 
(a gaseous sphere oscillating radially!), because, if ug 7 0 then this 
“sphere of gas would oscillate forever ($=0 71). Hence, this simple exam- 
ple shows that Stokes was incorrect on this point. Consequently, mea- 
Sured.values of secondary viscosity А turn out to be positive in direct 
-opposition to the Stokes hypothesis. Therefore, remember that for com- 
Pressible fluids in motion A» 0 and ув #0 


Viscous dissipation function then becomes $ = D: [(p-p)I + SI 


Фер: (p+terDl+sl=D: (pI +3] 


¢ = (ug Suc? «2 2:2 z0 


| 


Note that à= uy - چ‎ У 


Since 


ac 
"= 
" 
8 
* 
о 
+ 
Ө 
wn 
[n] 

— 
+ 
ю 
o 

=N 
N 
+ 
c 


23 


tr D 


ә 
" 
c 
* 
o 
* 
Ө 


and since 
2 2 


2 1 2 2 2 2 
Фуу + 072 + 033) = 30060) + (04705 + 05570," ED) 
we get the most general form for the viscous dissipation function (valid 
for all kinematically possible motions). Note that all the terms in © 


are positive! Hence, Ф 2 0 always! 


> > 


1 
of = ob + Ў[-р + (uç 20) 0] + 2 Ü (p) 


or 


Also 


а 


ONU 2 2 2 
Pasig esp) Schon Dy Dap e рур + Фуу") 


2,8 ес. 
аго Dig Dod] 


Ф = Au)? + zu (7a) + i6? - vii) 


26. 


Be aware that the usual values of "the coefficient of viscosity," и, 7 
appearing іп the common engineering literature, is actually a combined 
viscosity: tonto mud ا‎ 

vazu (= auth = 
This value came from the adopted linear Ps гө between the stress 


and the strain 


du 


Typus Bula 


If Ti are the stress tensor components (i,j = 1,2,3) and 
1 du, ч 
D.. = (к= Е LÀ are the components of the rate of deformation tensor 
ij D 8, x. 
(where u, are the velocity components), then 
T= [-p + (ug - 20) EDI + 2v 2 
becomes 
0, i#j 
5 PV pet İİ Dış em 
›і=ј 
For one-dimensional flows u, tu, = 0; uu; Til =T і 
so that A 
du 
T=-p+ (0 + 20)08 = prng a 


Hence, р = А +2. 
If the Stokes relation (A = - 2) holds, then combined viscosity is Б 
HL mm. 
e 
Both shear viscosity (и) and secondary viscosity (A) are actually 
not constants, but are functions of thermodynamic state. Ы 
Usually, the assumption that the temperature dependence of combined 


viscosity QU is given by a power law 


27. 


T 15 ; 
ue uo E where uy = const. 


is sufficient (although not exact!) for engineering computations. 

The exponent s has the value s = 0.647 for helium and s = 0.816 
for argon. These values both lie between the two most frequently used 
values provided by the kinetic theory of gases, namely, s = 0.5 and s 
1. 

Nevertheless, this simple-minded account of variability of the com- 
bined viscosity x is not acceptable if one wants to study shock wave 
structure where viscous dissipation and heat conduction are extremely 
important. The combined viscosity coefficient о) of undissociated air 


is well approximated by 


u = 1.462 a| چ‎ | 
с 


112 cm-sec — 
IX 


where T - absolute temperature in degrees Kelvin. 


The well-known formula of Lord Rayleigh suggests the value of exponent s 
to be s - 0.75 and claims its correctness for the temperatures between 


the freezing and boiling of water 


Me _ (T 40.75 
mim) 


o o 


Yet another formula of the same s-type is 


28. 


u l.5( T, + 223.2 
T 0 5 
a EN oe ОТЧ ES, 
Ho... (20 б 


уоп Міѕеѕ 


(Similar expression is valid for ©» 
t { v 


which;is valid in the limits between T = 32°F and T = 950°F.. 


u = 3.73 x 1077 [slug/ft.sec] = 0.000179 [g/cn.sec] 


at 15°C = 288°K = 518.4 °R 


But, in general, viscosity is the function of temperature, density and 
pressure, м 

сў He = ы. (Т, рур) + 
Therefore, one should include at least the effect of density on the 


combined viscosity 


4. du т 34 à Е z 
uo gu *a(o,T) = ї 


